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Q-i Abstract 



The discrete time quantum walk which is a quantum counterpart of random walk plays important roles in the theory 
of quantum information theory. In the present paper, we focus on discrete time quantum walks viewed as quantization 
of random walks on the path. We obtain a weak limit theorem for the time averaged distribution of our quantum 
walks. 



1 Introduction 

In the last decade, the discrete-time quantum walk (DTQW) has been extensively studied by many authors as 
a quantum counterpart of the random walk which plays important roles in various fields [Tll^Hl fTIIllTTll 1 4 II 2l?] . 
In the theory of quantum algorithm, quantum walk also plays important roles. In fact, quantum search 
algorithms are often reduced to quantum walks on the path, for example, search algorithms on the hypercube 
[TH] or glued binary tree [5] and an algorithm for element distinctness on the Johnson graph [3]. Therefore, 
investigations of DTQW on the path is beneficial. 

Rohde et al. [T7] studied periodic properties of entanglement for DTQW on the path determined by biased 
Hadamard coins numerically. Recently, Godsil [5] studied the time averaged distributions of continuous-time 
quantum walks (CTQW) on the path using the average mixing matrix. In the case of CTQW on the path 
determined by the graph Laplacian, the time averaged distribution consists of a uniform distribution and 
additional masses on the starting vertex and a vertex located at the symmetrical position about the middle 
point of the graph. In the present paper, we investigate the time averaged distribution of a DTQW. The 
quantum walk treated here can be viewed as a quantization of a random walk on the path. We should 
remark that the construction of our DTQW is closely related to that of Szegedy's work [19] . 

The rest of this paper is organized as follows. The definition of our DTQW is given in Sect. 2 and main 
result of this paper is stated in Sect. 3. The remaining section (Sect. 4) is devoted to the proof of our result. 
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2 Definition of the DTQW 



In this paper, we consider a discrete-time quantum walk on the path Pn+2 with n + 2 numbers of vertices 
Vn+2 = {0,1, ... ,n,n + 1} and the edges i?„+2 = {(*, « + 1) : i = 0, 1, ... , n}. For sake of simplicity of 
expressions, we choose Pn+2 because the number of the intermediate vertices in the; graph equals n. From 
now on we use n for suffix of operators although the number of the vertices of the graph equals n + 2. In order 
to define the DTQW, we use a Hilbert space Hn = span{|0, i?), |1, L), |1, i?), . . . , |n, L), |n, i?), |n + 1, L)} 
with \i,J) = \i) (El I J) [i G Vn+2, J G {L,R}) the tensor product of elements of two orthonormal bases 
: i £ Vn+2} for position of the walker and {\L) = -^[1, 0], \R) = ^[0, 1]} for the chirality which means the 
direction of the motion of the walker where "^A denotes the transpose of a matrix A. Then we consider the 
time evolution operator [/'^"^ on T-Ln defined by i7*^"^ = 5*^"^C^"^ with the coin operator C*^"^ and the shift 
operator S^") defined as follows: 



= |0){0| \R){R\ + li)01 ® (2|<^)(<^l - h) + \n + l)(n + 1| \L){L\, 



S'*")!^, J) = 



\i + l,L) if .J^R, 
\i-l,R) if J = L, 



where = ^\L) + ^^l-R) [Q < p,q < I, p + q=l) and J„ be the n x n idcntitiy matrix. 



Let X^' be the position of our quantum walker at time t. The probability that the walker with initial 
state is found at time t and the position x is defined by 



0/2) ([/("))* 1^) 



In this paper, we consider the DTQW starting from a vertex « G 14, and choose the initial chirality state 
with equal probability, i.e, we choose the initial state as |t/;)o = |0) ® \R) for i = 0, = \i) (g) \L) or 
= K) ^ \R) with probability 1/2 for 1 < i < n, or \'ip)n+i = \n + 1) §5 \L) for i = n + 1. For the sake of 
simplicity, we write Pi(Xj"^ = x) for P|^).(Xj"^ = x). We consider the time averaged distribution 



-(") (• \ 
Pi (x) 



lim E 

T-^oo 



T-1 



(n) 



X) 



t=0 



where the expectation takes for the choice of the initial chirality state. 



3 Main result 

Let be a random variable with distribution i.e., P(xj"^ = a;) = p^ix), we have the following 
weak limit theorem for X^ scaled by the system size n: 

Theorem 3.1 1. If the starting vertex i is fixed (finite) then 

-^(") 

^— ^Ci-5o + {l- a) ■ U{0, 1) (n ^ 00), 
n 

where U (0, 1) and Sq are the uniform distribution on [0, 1] and the delta measure at 0, respectively, and 

with the indicator function I a, i-e., Ia = 1> when A is true and I a = 0, otherwise. Here represents 
the weak convergence. 
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2. Let f{n) be a function such that rf{n) < n with < r < 1 anc 
vertex i is satisfing i/ f{n) — > r as n oo then 



f{n) —> cxD as n —> oo. // the starting 



C/(0,1) (n-^oo). 



The idea of the proof of Theorem 13.11 is based on [TO]. In this proof, the eigenspace of the following 
(n + 2) X (n + 2) finite Jacobi matrix Jn+2 induced by a random walk on P„+2 which is reflected with 
probability 1 at the boundaries and moves to the right and the left with probability p and q respectively, 
plays an important role: 



Jn+2 — 











V9 




pq 



pq 



pq 



/pq 



VP 







VP 





In addition, c,; in Theorem l3.1l is closely related to the stationary measure of the corresponding random walk. 
Let 7r„(i) (z = 0, . . . , n + 1) be the stationary measure of the corresponding random walk on Pn+2- By simple 
calculations, we have 



7r„(i) = 



2{l-{p/qj"^+'} 
2q{l~{p/q)'-+^}(P/l) 
2{T^W^^iP/lT 



if i = 0, 
if i = 1, . . . , n, 
if i = n + 1, 



for p ^ q and 



7r„(i) 




for p ~ q — 1/2. This suggests that co/2 = lim„_j.oo (0) and Ci = lim„_>oo 7rn(j) (1 < « < oo). In fact, 
7r„(i) is the same as square of the i-th component of the eigenvectors corresponding to the eigenvalue 1 and 
— 1 of Jn+2- These eigenvectors give the mass at the origin. 

It is well-known that the ballistic spreading of both CTQW and DTQW on the "infinite" line are 
characterized by weak limit theorems whose limit densities are quite similar except the extra term (1 — .t^)^^ 
in DTQW case [H]. We should remark that a relation between CTQW and DTQW on the infinite line is 
discussed in [15] . We also see more detailed discussion on the relations with respect to its weak limit theorems 
in [8] . In the case of CTQW determined by graph Laplacian on Pn+2 [h, , the time averaged distribution is 
almost the same as the uniform distribution. If we start from a vertex i e Vn+2, additional mass 1/ (2n + 2) 
is added to only two vertices i and n + 1 — i. Because the additional masses tend to zero as n — ?■ oo, the 
limit distribution of CTQW corresponding to Theorem 13.11 is the uniform distribution U{0, 1). Therefore, 
the localization at the origin in the sense of Theorem 13.11 is not observed in the case of the CTQW. On the 
other hand, we can see in [7] a weak limit theorem for a DTQW Xj""-* on the half line starting from the 
origin whose system is infinite in advance. In this case. 



X. 



(oo) 



Co • (5o + (1 - Co) • K{2^) {t oo) 



where K{r) (0 < r < 1) has the following density: 

x^ 

f{x;r) 



1-VT 



fK{x;r) ■ I[o,oo]{x), 



where /if (cc; r) is the Konno density function [T^I13) . This suggests that infinite nature of underlying graph 
changes the shape of limit distribution from a uniform distribution. It is a future problem finding connection 
between CTQW and DTQW on finite graphs in the sense of [15^ . 
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4 Proof of Theorem 13.1 



In order to prove Theorem 13. 11 we consider the eigenspace of the Jacobi matrix Jn+2 at first. 
Lemma 4.1 1. Let A™ (to = 0, 1, . . . , n, n + 1) be the eigenvalues of the matrix Jn+2- Then we have 

Ao = 1, 

Xm = 2^^008 6*™ (to = 1, . . . ,ri), 



-1. 



where 9m = nnr/{n +1). 

2. Let Vm be the eigenvector corresponding to an eigenvalue Xm for each < to < n + 1 . Then we have 



T ~ 
vo = 



p- q 



1 

1,— , 



(v^) (v^) 



Vq \/9 



— Cm ■ 



1, 



Xm. I XmUi 



(jri) 



where C„ 



Vq \ Vp 

1 



' V9 

Vq Vp 



-,(vW^) 



(7 - = and cos(y9m — Xm, for 1 < to < n, 

sm 9m. 



2p sin^ 



T ~ 
V„+l 



p-q 



1, ^, ^ ^ ^ , . . . , — , (-1)"+^ [VWq] 



Vq Vq 



The eigenvalues 1 = Aq and — 1 = A„+i are distinguished from the other eigenvalues Xj (j = 1, . . . , n), that 
is, the eigenvalues Xj (j = 1, . . . , n) are continuously distributed while the eigenvalues 1 and —1 remain as 
mass points when n — > oo. This structure of the eigenspace affects the weak limit theorem. 

Proof of Lemma 14.11 

Although it is well-known result (see e.g. but we show the proof because of improvement of the 
readability. Let Dk — det(A/fc — Jfc) and Ek = det(A/fc — VPq^k) — iVPq)'^ det(A/y^/fc — Ak), where Ak is 
the adjacency matrix of the path Pk- We can easily obtain the following recurrence equation: 



Dk+2 = A~£^fe - XEk-i +pqEk-2- 
On the other hand, we have the following recurrence equation by a simple calculation: 

Eo/iVpqf = i> 
EiKVpqf = >^/Vpq, 
[x/Vm) Ek/iVrnf = Eu+iKVmf^' + Ek-i/iVm)"-' (fc = 2, 3, . . .). 

Noting that |A| < 1 by the Perron- Frobenius theorem since < p,q < 1 and p + q = 1, this implies 



(4.1) 



Ek = 



'a'=(/+i -g'=+i)/(?j-(7) if|A| = l, 

iVpq)''Ukix/Vm) if |A| < 1, 



(4.2) 



where Uk{x) is the (monic) Chebyshev polynomial of the second kind, i.e., Uk{x) = sin(fc + 1)9/ sm9 with 
a; = 2cos^. Combining Eq. (|4.ip with Eq. (14.21) . we have Dk+2 = for |A| = 1. For |A| < 1, we have 



Dn+2 = iVmnx' - 1) ''"^"V^^ ^ith = A/(2^). 

sm 9 
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In this case, the solutions of Dn+2 — are 9 — rmr/ {n + 1) (m = 1, . . . , n). Therefore we obtain the desired 
result. 

Next we estimate the eigenvector. Let Po{X) — 1, Pi{X) — A, -P2(A) = det(A/2 — ^^^2) XPi{X) — qPo{X), 



Pfc(A) = dct 







pq 







pq 



fq X 



{2<k<n+l), 



where the size of the matrix defining Pk{X) is k x k. Then we have the eigenvector v„i ~ v„i/||v,„|| (see 
Lemma 1.91 of where 



^v,„= [Po(A„0,A(A™)/Vg,^2(A„0/V9(;^,---,^n+i(Am)/V9(MF ■ 
For 3 < fc < n + 1, we have Pk{Xm) = XmEk-i — <lEk-2{Xm)- Combining this with Eq. (I4.2p . we obtain 



" if TO = 0, 

Pk{Xm) = ^ A„(^)'=-l[/fc_i {Xm/^/pq) - q{y^)''^^Uk-2 {Xm/y/Pq) if 1 < TO < n, 

_^p-^k-i if TO = n + 1. 



Finally, we calculate ||vm|p. By a direct calculation, we have 



p-q 



Let f/^™'' = Uk{Xm/ \/pq) — sin(fc + l)0„j/sin0m with dm — rmr/{n + 1) for 1 < to < and < k < n. 
Noting that u'^^ = ut"^ = 0, we obtain 



p>L + g 

pq 



■n—i 

E (^r ) 



\/g Vp ) 



2 0\2 "^1 



J=0 



Using (A„/v/pg)C/j™^ = L/]!;^) + C/j™^ and ui"l\ - (-1)™, we have 



. 2 n— 1 

v„,|2 = l + ^ 



9 -2 

^ sm 93 



3=0 



3=0 



with sin(/3„i = Xm- On the other hand. 



■ 2/ ■ 



sin^(j + i)e„, _ 1 



3=0 



3=0 



sm 



2 sin^ 



n - ^ cos 2j9m ^ = 



2 sin^ 6*. 



^In-fy^") cosn0,4 



2 sin^ 9r, 



This completes the proof. □ 
Next we obtain the eigenspace of the time evolution operator C/'"^. Note that there are 2n + 2 numbers 
of eigenvalues by the definition of the Hilbert space defining our DTQW. 
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Lemma 4.2 1. Let fik {k — 0, ±1, ±2, . . . , ±n, n + 1) be the eigenvalues of the time evolution operator 
[/("). Then we have 

M±m = 6='=*'^'" (m = 1, 2, . . . , n), 

where i = \/— T and cos if „i = Kn G Spec (Jn+2) \ {I7 ^1} with (p,,,, G (0,7r). 

2. Let Ufc be the eigenvector corresponding to an eigenvalue fik for each k = 0, ±1, ±2, . . . , ±n, n + 1. 
Then we have 



n+l 

3=0 



with 



= j-*-(c/f\-/i..y^(7j!y »/l<J<" + l, fc = ±l,±2,...,±n, 

otherwise, 



(p/g)" + i-2g 



[t^k Vp/Qy ^ j/ 1 < j < n + 1, fc = 0, n + 1, 




^^^^1 {p/qr+'-2q • (MfcV-P/g) ifO<j<n, fc = 0,n + l, 



= i • f/ifet/f ^ - v/g/p ) / \/ ? < J < n, fc = ±1, ±2, . . . , ±n, 

ot/ieraise. 

The eigenvalues of the Jacobi matrix J„+2 are projections of the eigenvalues of C/*^") to the real axis from 
the unit circle of the complex plane. Therefore the eigenvalues (of f/*^"^) 1 = /io and — 1 = /i„+i are also 
distinguished from other eigenvalues (of C/^")) (j = 1, . . . ,n) because the eigenvalues (of Jn+2) 1 = Aq 
and —1 = Xn+i are distinguished from the other eigenvalues (of Jn+2) Xj (j = 1, • ■ • The eigenvectors 
corresponding to the eigenvalues 1 and —1 give the mass at the origin in Theorem 13. II 

Proof of Lemma 14.21 

Let Vm(j) (j — 0,1, ... , n, n + 1) be the j-th component of the eigenvector corresponding to Am € 
Spec (J„+2). We define the following vectors: 

n 

a™ = v„(0)|0, R)+Y, v,„(j)|j) ® (Vg|i> + v^|i?)) + v„(n + l)|n + 1, L), 

It is easy to see that [/^"•'am = b,„ and [/("'b™ — 2A„jb„j — a,„. This implies that C/^"^ is a linear 
transformation on the Hilbert space spanned by the two vectors a^ and b^. Moreover, the inner product 
(am,bm) equals Am||am|| ||bm|| = Xm- Therefore, for m = (resp. m = n+l), i.e.. A™ = 1 (resp. A™ = -1), 
case, we have b^ = Amam- Thus Aq = 1 and A„+i = — 1 are eigenvalues of U^^^ with corresponding 
eigenvectors ao and a„+i, respectively. For m ^ 0,n + 1 (Am ^ ±1) case, {am,bm} is a base of the Hilbert 
space spanned by am and bm- The representation matrix of C/^"^ with respect to the base {am,bm} is 
[1 2Ai]- Therefore, the eigenvalues are e^*"^™ and the corresponding eigenvectors are a™ — e^^'^'^bm- It is 
easy to see that ||am, — e^"^™bm||^ = 2(1 — A^) = 2 sin^ ipm. Combining these facts and Lemma HTTl we have 
the explicit form of the eigenvectors after a long but simple calculation. □ 
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Now we estimate the distribution pf"^ of the random variable X j""* . By the assumption of the choice of 
the initial state, we have 



lim 



T-1 

limT^oo ^ Y^^Q 



lim 



T Z^t=0 



((x|®/2) (|0)®|i?)) 

5Ej=l,h||((^I®^2)(C/("))'(N>®|J)) 
(8/2) (|n + l)(»|L)) 



if i = 0, 

if 1 < i < n, 

if i = ?i + 1. 



Using the spectral decomposition = ^|,Ufeuj^ and limT->oo(l/r) Y^=o s'^* — <5o(^) (mod 27r), we 

obtain 



if i = 0, 

if 1 < i < n, 

if i = n + 1 , 



because all eigenvalues of J7^"^ are nondegenerate. Note that it is enough to treat the first two cases of 
p["\x) in our setting. From Lemma [4.21 we have the following expressions: 



Po U) = 



So{j) + {l-So{j))-{p/qy-' 



+ 



2p 



E 



sin 6„ 
sin 



p (C7j")) ' - cos 2(p™ (tj^Ziy + q (U^-!^ ' 



(4.3) 



for i = and 



{p - qf 



q{{plqr+^-2qY 



{p/qr-H5,{j) + {l-6,{j))-{p/q) 



J-1 



smf 



sm Lp„ 



P 



,s2^„(t/]"))%,({/ji)' 



for 1 < i < n. 



Prom now on, we estimate P(xj"'' < an) = J2j=oPi"^'ij) (0 ^ 1) where [zj denotes the largest 
integer not greater than z. At first, we consider the first terms in Eqs. (|4.3p and (j4.4p . We can easily see 
that 



COS 



(4.4) 



E UoU) + (1 - <5o(j-)) i(p/<z)-^| = 1 + ^ X ^ ^""^'^ 



[anj 



ip/q) ■ 

Therefore the first terms converge to zero ii p > q and finite i case and infinite i case. For p < q and finite i 
case, the first term in Eq. (j4.3p converges to 1—p/q, also that of Eq. (j4.4p converges to {l—p/q){p/qY^^ / (2q). 
Next we estimate the second term in Eq. (|4.3p . In this case we should consider the following summation: 



[anJ 



^2 E E 



sin 6„ 

sin if,- 



p {u^"'^y - cos2^,„ {u^:^iy + q (u^ziy 



j—O m — 1 

Note that this summation is finite thus we can exchange the order. Using the following relation: 

I 



(m) 



2 sin 



^ |; - cos(?0„){//™n , 
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we obtain 



m—1 j—0 ^ ^ ^ 

^ (S^) ~ cos(2^™)} { lan\ - cos(LanJ0„)i/f"jj } 

m—1 



m—1 



sin 9m \ ^ 1 



sm iy9m / n 
sin fc \^ dk 



\ sin / TT ' 

where f{n) ~ ^('t^) denotes lini„_>.oo f{n)/g{n) — 1 and cosiy9(fc) = 2^/pgcosfc. This imphes the second term 
in Eq. (|4.3p converges to the uniform distribution C/(0, 1) with total mass 2p (sinfc/ sin(p(fc))^ dfc/7r. This 
total mass can be calculated directly but it is obvious that this integral equals 1 — cq. 
We can estimate the second term in Eq. (j4.4p by using similar argument as follows: 

m—1 ^ / 



X [p (t/i™^) ' - COS + 9 {my 

^ psin{(i 4- l)fc} — cos2i^(fc) sin(ifc) + qsin{(i — l)fc} dk 
sin(^(fc)2 TT 

Generally, it is hard to calculate the explicit value of the last integral but it is obvious that this integral 
equals 1 — . This completes the proof. 
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